Central mode and spin confinement near the boundary of the superconducting phase 

in YBa 2 Cu 3 6 .353 (T c =18 K) 
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We have mapped the neutron scattering spin spectrum at low-energies in YBa2Cu3 06.353 (T c =18 
K) where the doping ~0.06 is near the critical value (p c =0.055) for superconductivity. No coexistence 
with long range ordered antiferromagnetism is found. The spins fluctuate on two energy scales, one a 
damped spin response with a ~2 meV relaxation rate and the other a central mode with a relaxation 
rate that slows to less than 0.08 meV below T c . The spectrum mirrors that of a soft mode driving a 
central mode. Extremely short correlation lengths, 42±5 A in-plane and 8±2 A along the c direction, 
and isotropic spin orientations for the central mode indicate that the correlations are subcritical with 
respect to any second order transition to Neel order. The dynamics follows a model where damped 
spin fluctuations are coupled to the slow fluctuations of regions with correlations shortened by the 
hole doping. 

PACS numbers: 74.72.-h, 75.25.+Z, 75.40.Gb 



Compared to the extensive studies of spin fluctua- 
tions deep within the superconducting (SC) phase of 
cuprates^^ less is known for low doping. There is little 
consensus as to the nature of the precursor phase from 
which the SC phase emerges when holes are doped into 
a planar array of correlated S=l/2 Cu spins. Whereas 
superconductivity in I^-zSr^CuC^ (LSCO) emerges 
from a non-metallic phase, the precursor phase for 
YBa2Cu306+a; (YBCOg+a;) remains controversial 3 - 4 - 5 - 6 . 
Probing the spin structure in the vicinity of the bound- 
aries of the antiferromagnetic (AF) and SC phases is nec- 
essary to determine whether antiferromagnetism and su- 
perconductivity compete, coexist, or are separated by a 
novel phase. 

For single layer LSCO, an insulating phase with frozen 
short-range magnetic order separates the AF ordered 
from the SC phaseS&SiiS. At the SC boundary, the mag- 
netic structure abruptly changes by a 45° rotation of the 
incommensurate wave vector of static spin stripes, so the 
transition to superconductivity is first order. In contrast, 
the temperature dependence of the commensurate scat- 
tering in La2Cu04 lightly doped with Li is interpreted 
as indicating the proximity of a quantum critical point 
(QCP)ii. It is clearly important to study other cuprates 
such as YBCOg+a; to determine the nature of the state 
near the edge of the SC phase. 

For the bilayer system YHGOq+ x oxygen doping oc- 
curs in the chains well away from the planes and complex 
structural changes as in LSCO are not evident. There 
have been suggestions that the AF and SC phases may 
meet at a pointiSii^. Recently Sanna et al~ claimed 



from fj,SK data that the ordered AF state coexisted with 
superconductivity with a Neel temperature exceeding 
that of a freezing transition to a 'cluster spin-glass' state. 
This contrasts with LSCO and Yi_j / Ca 2/ Ba2Cu306+ a ; 
where AF and SC phases are well separated as shown 
in Fig. 4 of Sannaii 5 - 

Another unique aspect of YBCOe+z is that the SC 
phase reveals a resonance peak whose energy scales with 
T c for at least a limited range of dopin g 4 - 2 d 6 . Quantum 
critical theories suggest that the resonance peak may be 
a soft mode of the SC phase— and therefore it is crucial 
to study spin behavior close to the critical doping for 
superconductivity of p c =0.055^ to test whether the anti- 
ferromagnetic QCP lies within the SC phase or whether 
it is contiguous. 

Here we approach the SC boundary by studying 
YBCO6.353 where T c =18 K is suppressed to one fifth of 
optimal. The hole density of p^0.06 exceeds the critical 
value (p=0.055) by only ~9%. Earlier results^ 8 - with a 
similar stated oxygen concentration gave a much higher 
T c and so lie much farther from the critical doping. We 
will show that just beyond the critical doping the spins 
fluctuate on two distinct time scales, one giving a broad 
energy feature and the other a central mode at zero en- 
ergy. A single spectral form is found in which damped 
spin fluctuations transfer strength to an intense central 
mode with an extremely narrow energy width similar to 
the spectrum near antiferroelectric transitions 1 ^ .. The 
absence of long-range spin order and short spatial cor- 
relations indicate that the system may lie farther from 
a second order phase transition to long-range AF order 
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FIG. 1: (Color online) Polarized spin-flip neutron scatter- 
ing, a) The magnetic spectrum along [1,1,0] measured by 
subtracting the VF from the HF data. Both inelastic 6) and 
elastic scans c) show the spin directions are isotropic since 
the HF intensity is twice that for VF. The horizontal bar in 
c represents the resolution full- width. 



FIG. 2: (Color online) Elastic scans a) Comparison of short- 
range incipient 3D correlations normal to the planes and the 
coupled bilayer model (line). An 80 K background has been 
subtracted, b) and c) Polarized neutron spin-flip scattering 
for horizontal and vertical field. The fact that HF=2VF in- 
dependent of L shows that the scattering is isotropic. The 
enhancement of the L=l peak for E/ = 14.6 meV relative to 5 
meV arises from resolution effects. 



than the lower critical concentration for the SC phase. 
Taken with the conservation of spin orientational sym- 
metry found with polarized neutrons, the spins display, 
not a freezing transition to a cluster phase, but a grad- 
ual formation of only short-range correlations that are 
intrinsic to the SC phase. 

Thermal neutron measurements were made with the 
DUALSPEC spectrometer at Chalk River Laboratories 
using PG(002) as monochromator and analyzer and a PG 
filter in the scattered beam. The horizontal collimation 
was set to [33' 29' S 51' 120'] and [33' 48' S 51' 120'] for 
E/=14.6 meV and 5 meV, respectively. The setup for po- 
larized beam measurements was described previously^ 
Cold-neutron measurements were made with the SPINS 
spectrometer at the NCNR with a PG (002) monochro- 
mator and analyzer. For E/=2.9 meV the collimation 
was [guide 80' S 80' 80']. For energy transfer less than 
2 meV a Be filter was placed before the monochromator 
and above 2 meV a BeO filter was placed in the scattered 
beam. 

Seven ^1 cc crystals were co-aligned in the (H,H,L) 
scattering plane with an overall mosaic of ^1.5°. From 
DC magnetization (Fig. 3c) the critical temperature was 
measured to be T c =18 K with a width of ^2 K. The 
orthorhombic lattice constants were a=3.843, 6=3.871, 
and c=11.788 A. Superlattice peaks associated with CuO 
chain ordering were not observable indicating that oxy- 
gen ordering in the chains is short-ranged. 

Polarized neutron measurements of the HF-VF sub- 
tracted spin-flip scattering (Fig. ^,) probe the spin fluc- 
tuation spectrum along the [1,1,0] direction. They show 
that the spins fluctuate on two energy scales, one showing 



a resolution limited central peak (Huj=0) with very slow 
dynamics, and another with a faster timescale giving a 
broad inelastic peak at ^2 meV (see also Fig. 0J. 

Scans through the AF wave vector along [H,H,2] (Fig. 
HJd, c) and along [|,|,L] (Fig. c) show that the HF 
intensity is twice the VF intensity. It follows that the 
fluctuations are equally strong within and normal to the 
(H,H,L) scattering plane, even as the varying L changes 
the amount of [0,0,1] fluctuations sampled. This shows 
that the spins are isotropic in direction for the near-static 
central mode. Fig. ^ shows that inelastic fluctuations 
for Q=(|,|,2) are also isotropic. At all temperatures 
the spin rotational symmetry is conserved. This state 
clearly differs from the anisotropic critical fluctuations 
that would have been seen were the system approaching 
an AF phase, in which moments aligned in-plane would 
have caused an energy gap for out-of-plane fluctuations 
polarized along [001]. The low-temperature isotropic 
magnetic ground state observed in YBCO6.353 is clearly 
dramatically different from that of the insulator. The 
static spin cross-section as shown later is dynamically 
coupled to spin excitations and so is not the result of a 
phase-segregation of holes or a doping gradient. If this 
had been so, the clusters would have mimicked the prop- 
erties and anisotropy of the antiferromagnetic insulator. 
It is presently not clear how this ground state relates 
to that suggested to explain NMR results in LaCu04+<5, 
crystal field excitations in ErBa2Cu30;c, and /iSR data 
in LSCO and Ca doped YBCOi 15 i 21 i 22 

From scans with E/=2.9 meV along [110] we extract 
an in-plane correlation length of 42±5 A for the central 
mode from a resolution-convolved fit to a lorentzian. 
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The integer-centered peaks along L in Fig. [2^ reflect 
the lattice periodicity, but have a large width. The elastic 
scattering is well described by a model where internally 
AF bilayers are weakly coupled along [001]. The struc- 
ture factor along L for fixed in-plane wave vector (5)0) 
is proportional to 
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with £ the correlation length and d the bilayer separa- 
tion (d/e=0.28), and to the anisotropic Cu 2+ form fac- 
tor—. The resolution convolved fit (solid line) gives a dy- 
namic correlation length of only 8±2 A along [001]. This 
interplanar correlation of less than one unit cell shows 
that the spins are primarily correlated within the two- 
dimensional (2D) CuC>2 bilayers. The short-range na- 
ture of the spin correlations indicate that YBCO6.353 is 
strongly sub-critical with respect to 3D long range AF 
order even at low temperature. The presence of short- 
ranged elastic correlations and two distinct timescales 
differs substantially from the behavior near a second- 
order phase transition where the energy lincwidth is di- 
rectly related to and scales with the correlation length. 
We emphasize that short-ranged spin correlations do not 
imply the formation of physical clusters except in the 
most naive interpretation; a uniform state exhibiting su- 
perconductivity and finite spin correlations is a possible 
ground state. 

Fig. |21 displays the peak intensity and energy width of 
the central mode spectrum from Gaussian fits. The peak 
intensity grows on cooling, dramatically and smoothly 
(panel a) , with no sign of a transition, and is unperturbed 
by the sharp SC transition at 18 K (panel c). On cooling, 
the central mode energy width declines to the resolution 
limit of ^0.08 meV (panel b). A possible picture is of 
spins coupled through a fermionic response associated 
with holes. 

The appearance of a sharp central mode, whose energy 
scale is much smaller than that of a damped inelastic 
mode, bears a strong similarity to the soft and central 
mode spectrum in SrTiC>3. Halperin and Varma— and 
Shapiro et aZ— showed that a coupling of a damped har- 
monic mode to a slowly fluctuating object would give 
rise to a central mode. We surmise that the object in 
YBCO6.353 results from the slow fluctuations of the many 
spins within a correlation range frustrated by the metallic 
hole doping. 

The spectrum at (tt, it) (Fig. 0} was obtained from the 
difference between constant-Q scans at the peak of the 
magnetic scattering, Q=(i,i,2), and at a background 
position (0.35, 0.35, 2). The solid lines are the result of 
fits to the phenomenological modeliS^ 
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FIG. 3: (Color online) Panel a) illustrates the growth on cool- 
ing of the peak intensity of the central mode from a Gaussian 
fit and compares it with the volume fraction of static spins 
measured by Sanna et al. using muons. Panel b) shows how 
the full-width of the central peak drops smoothly to the res- 
olution limit (broken line) below T c (dotted line) and panel 
c) shows the sharp Meissner transition at 18 K. 



which describes a mode of renormalized frequency fin and 
damping T coupled to an object that is hopping slowly 
at a rate 7=1/7-. The Bose factor is [n(uo) + 1] and 8 
is a measure of the coupling between slow and fast spin 
modes. The rapid growth of the central mode intensity 
on cooling (Fig. |3Ji) then arises, for fixed ?ir=10 meV 
and amplitude, from a softening of the harmonic mode 
On and from a slowing of the hopping rate r _1 . The 
model of a soft mode driving a central-mode gives a good 
account of the two energy scales in Fig. 0]and describes 
the data over a broad range of temperature and energy. 

In addition to the central mode, which is absent in 
YBC06.3^, the spectrum of Figs. ^ and 0] shows that 
the spin energy scale, defined by the ~2 meV peak in 
S(Q,uj) 1 has decreased by an order of magnitude in 
YBCO6.353 (Fig. ^) from the 33 meV resonance en- 
Qj-gyi^M Q £ YBCO6.5. This is larger than the doping 
ratio 0.06:0.09 and may indicate the predicted collapse 
of the spin energy on approach to the SC boundary— 

In recent /^SR measurements, Sanna et al*^ inferred 
that a freezing transition at T/ to a 'cluster spin-glass' 
state occurred within the SC phase. A direct compari- 
son with the neutron results (see Fig. [2J shows different 
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YBC0 635 , E f =2.9 meV 
Q=(1/2, 1/2. 2.0) 

T=2 K 




FIG. 4: (Color online) Temperature evolution of the back- 
ground subtracted spin spectrum at Q=(7r,7r). The line shows 
the response of a damped mode that on cooling softens and 
drives up the intensity of a central peak arising from slow 
hopping of a spin cluster (see text). The parameters for the 
fits at T=2, 20, 30 K respectively are M7 =0.9 ± 0.3, 2.4 ± 
1.1, 3.4 ± 1.9 meV; r/h= 14.9 ± 3.0, 9.2 ± 1.2, 1.8 ± 0.8 
meV" 1 ; ft6= 3.9 ±1.3, 2.0 ±0.9, 2.0 ±0.9 meV. 



onset temperatures for the different probes. This illus- 
trates that T f (measured with muons) represents a slow- 
ing down beyond the longest muon timescale. The freez- 
ing temperature does not represent a sharp transition. 
The muon phase diagram^ also suggests a coexistence 
of the SC phase and Neel order. Our results prove there 
is no coexistence of long-range AF spin order just inside 
the SC dome. Instead neutrons show, directly, that spin 
correlations are short-ranged and isotropic and that there 
is no phase transition to a 'glass' or 'cluster' state. 

We believe the slowing down of /uSR in YBCO arises 
from the same central peak as observed here in supercon- 
ducting low- doped YBCO. 

The different temperatures where the slow dynamics 
becomes observable (Fig [3J reflect the timescales of the 
probes. These results do not require the existence of 
literal clusters or of phase separation of holes. Unlike 
neutrons, /iSR cannot directly detect both spatial and 
temporal information of the spin correlations, but com- 
parison shows that the muon data is consistent with the 
presence of the two energy scales and short-ranged cor- 
relations observed here. The short-range correlations co- 
exist with the nascent superconducting ground state of 
YBCO - a complex ground state with slow fluid-like dy- 
namics 

We have shown that in the superconducting state of 
YBCO a central mode arises from fluctuations that grow 
dramatically and become extremely slow in the low- 
temperature limit in the presence of faster scale of relax- 
ing spin excitations. The short-range correlations within 
and between planes show the ground state consists of sub- 
critical 3D-enhanced spin correlations. The isotropic spin 
orientations suggest that the spin configuration is radi- 
cally changed from the antiferromagnet by the relatively 
small hole doping. 
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